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Some properties of an ideal gas of massive bosons in a mean field potential 
and, confined between two infinite parallel slabs in a d-dimensional confignra- 
tion space are investigated systematically. Here, one particle density of states 
approach is employed to stndy the critical temperatnre, shift of density, Casimir 
effects and critical exponents, starting from the evalnation of the grand canoni¬ 
cal free energy in d-dimension. We have fonnd that, the shift of density, Casimir 
force and the critical temperatnre depend on the space dimensionality. Bnt the 
Casimir force decays as an inverse power law of the distance between two slabs 
in the condensate and, decays exponentially in the non-condensed state sitnated 
very close to the point of phase transition. Most importently, this stndy enabled 
ns to predict the shift of the density of excited bosons dne to mean held poten¬ 
tial, and also the dimensional dependence of the critical exponent. The form of 
the critical exponent is fonnd to be for the imperfect Bose gas. This leads 
to a valne of critical exponent 1 for d = 3. 
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I. INTRODUCTION 


After the demonstration by Einstein that there is a possibility of condensation ( 
which is referred to as the Bose-Einstein condensation) of the ideal Bose gas (IBG) many 
attempts have been made to understand this phenomenon. Since then the bulk behavior 
including condensation characterized by the thermodynamic properties, occupation of 
the states and their fluctuations are investigated[1-5]. But after the observation of the 
Bose-Einstein Condensation (BEC) of cold atoms in a trap in 1995 [6-8] the interest of 
better understanding of Bose condensation has been strongly stimulated. As a result 
different studies are made on a Bose system analyzing the effects of potentials [9-14], 
interaction between particles [15-24], space dimensionality [25-30], criticality near the 
phase transition [26,31,32], and hnite particle number [9]. 

Recently, some attempts have been spent to study imperfect Bose gas (IMBG) plug¬ 
ging in the mean held potential[15,33-35]. Napiorkowski and Piasecki, in [33], evaluated 
the grand canonical potential considering the 3-dimensional space along with the mean 
held potential, within the perspective of the Quantum Statistical Mechanics. However, 
the density of states for this type of imperfect Bose gas has not been studied yet, al¬ 
though it plays an important role in statistical physics [14,25,27]. On the other hand, 
theoretical studies demonstrate that the space dimensionality and the type of poten¬ 
tial ahects diherent properties such as the density [20], critical temperature [20], the 
thermodynamic Gasimir force and the critical exponents [33]. It is, therefore, inter¬ 
esting to derive one particle density of states in d-dimensional conhguration space for 
the imperfect Bose gas (IMBG) that involves the mean held potential. It is also very 
much interesting to examine then if this density of states approach can predict different 
physical properties specihcally the Bose condensation, density shift and Gasimir effects 
of the IMBG correctly. 

In a real system, interaction among the particles does vitally exist, but it makes the 
problem difficult to solve analytically. In order to make the problem solvable, and to 
retain the essential physics intact one can represent a real system approximately by an 
ideal system of non interacting particles in the presence of a potential or a trap. The 
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validity of the approach is nicely justihed by the experimental fact that the influence of 
interaction between particles on the Bose-Einstein condensation transition temperature 
is only of several percent [16]. 

In this report, we intend to derive hrst the one particle density of states for the 
ideal Bose gas in a mean held potential (MFP) using a d-dimensional conhguration 
space and then to apply it to evaluate the grand potential energy for a system conhned 
between two parallel inhnite slabs separated by a distance of D. The magnitude of D 
must be much much larger than that of the thermal wave length of the massive bosons. 
Secondly, the shift of density in the excited states of the imperfect Bose gas due to the 
MFP, and some other characteristic properties across the phase separating point will 
be studied. Finally different features of the Casimir force and critical behavior will be 
examined with some rigorousness. We should mention here that, the reliability of the 
derived analytic expressions describing different properties of the imperfect Bose gas 
(IMBG) should be justihed by comparing with other known results. To this end, we 
intend to examine the limiting cases as the MFP tends to zero (i.e. a —)■ 0). We note 
here that, in the limit a —)■ 0 all results (presented below) show complete coincidence 
with those of the ideal Bose gas as expected. 

The magnitude of the Casimir effect depends on the boundary conditions (be) used. 
So, we intend to use the Neumann boundary condition explicitly in the theory and then 
to discuss the possible impact of other boundary conditions on the results. We note 
here that the sums for the periodic, Neumann and Dirichlet boundary conditions starts 
at — cxD, 0 and 1, respectively. However, the latter two bes can be expressed in terms 
of sums from —oo to -|-cxd. On the other hand, the operator relations for Neumann 
and Dirichlet bes are —)■ respectively. The 

periodic be is idealistic in the sense that the corresponding surface energy density 
becomes zero, while for others it becomes nonzero. We also note that the magnitude of 
Casimir interaction energy is the largest for periodic be and smallest for Dirichlet one 
| 22 ]. 

This report is organized in the following way. Derivation of the one particle density 
of states in d-dimensional space, grand canonical free energy, density shift, Casimir 
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amplitude and force are made in section 2. The critical behavior is also examined in 
the same section. Section 3 is devoted to present the results and discussions. This 
report is concluded in section 4. 


II. THEORETICAL EVALUATION 

A. One particle density of states in d-dimensional configuration space 


The repulsive pair interaction between a pair of identical massive bosons can be 
described within the mean field theory by a/V where a is a positive constant and V 
is the volume of the Bose gas. Let us now look at a particular boson in the N boson 
system which is moving in the mean held due to the rest of (N — 1) particles. Obviously, 
the average potential energy experienced by the tagged boson is y(N— 1) ^ yN. The 
one boson Hamiltonian will, therefore, be 

2 

H = ——h an (1) 


where p is the momentum of the boson and m is its mass; n denotes the number density 
of the bosons {N/V). 

Let us assume that, the bosons are enclosed in a d-dimensional volume (denoted by 
U('^) = L being the edge of the rectangular box and L —)■ oo. The spacing 
between energy levels will therefore be very small, so the summation over the states 
can be replaced approximately by an integration. Therefore, the bulk density of states 
in the phase space with spatial d-dimension is 
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e in equation (2) is the energy eigen value of the Hamiltonian operator of equation (1). 

We now assume that the Bose gas in the mean held potential (henceforth will be 
referred to as the imperfect Bose gas, IMBG) is placed within two inhnite parallel 
slabs such that the d-dimensional volume D. The density of states in the 
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((i-l)-dimensional surface is now 


7i(e) 


y(d-^) 


m 1 

' r(^^) 



B. The grand potential of the system 


(3) 


As D is finite, the spacing of the energy levels will be large along the dth direction. 
So, the summation over energy levels along the dth direction cannot be approximated 

by integration. Then the grand potential (with wave number kd = ^ 1; I = 0,1, 2, 3.), 

may be expressed applying the Neumann boundary condition as 
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where = exp(/5/e),/5 is inverse temperature times Boltzmann constant, r(4^, —(danr) 
the lower incomplete gamma function and A = hjy^2wmk^, h being the Planck’s 
constant. Now separating the I = 0 term and using the following identity relation 
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here, we have assumed that 
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Equation (7) is valid when the magnitude of lower limit of the integration in the incom¬ 
plete gamma function is very small; in other words, when the strength of interaction 
is weak. We note that the lower limit contains a sum over r. So a question remains 
whether the approximation is valid for large value of r or not. The answer is yes be¬ 
cause from equation (4) we see that the lower limit of the incomplete gamma function, 
—/danr becomes large for large r and changes the value of gamma function to some 
extent. But, at the same time the multiplying factor in the term goes to be 

very small relatively. As a result, the contribution of each term to the energy at large 
r is insignihcant despite the ratio of gamma functions in Eq.(7) deviates from 1. How¬ 
ever, for the present study, equation (6) is the hnal expression for the grand potential 
that involves the d-dimensional conhguration space and the mean held potential (MFP) 
energy. 

In the absence of the MFP energy i.e. with a —)■ 0 equation (6) stands as 
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where ga{,z) is the Bose-Einstein function [36,37]. For d = 3, equation (8) becomes as 
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In equation (9) the bulk energy density is —ksT gs (z) /A^ and the surface energy density 
is —g 2 {z)/\^ when fj, < 0 and, these densities in the condensate (/i = 0) become 
-ksTCiD/X^ and -k bT({2), respectively. The third term on the right hand side is 
known as the Casimir interaction term 
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where duiD = If is noticed that, equation (10) completely agrees with the 

properties of ideal Bose gas (IBG) [38,39]. This result, however, conhrms that, the 
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first and the second terms in the d-dimensional expression ( Equation (6)) are surface 
and bulk energy densities, respectively. The third term is nothing but the Casimir 
interaction in the d-dimensional conhguration space for the IMBG. 


C. Critical Temperature 


The second term on the right hand side of equation (6) describes the bulk energy 
density 

ksT 
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■ Q d+2 (^Z ) 


( 11 ) 
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where z' = From this knowledge of the bulk energy density in d-dimension 

one can easily evaluate the bosonic number density in the bulk, 

dub 
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So, one can write 


(ksT)U(l: 


(13) 


where ( denotes the Riemann zeta function. Now at the critical density one can write 
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which leads to the critical temperature in d-dimensional space 
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From equation (15) one can easily show that the fraction of Bose particles in the ground 
state, 
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where n is the sum of the number of bosons per unit volume in the ground and excited 
states i.e. n = no + 

This is of great interest to see the affect of the MFP on the shift of density of the 
excited bosons. In order to examine it we consider the lowest order of MFP to have 
z' = e^^(l — I3an). For the MFP to be weak I3an « 1, so it is obvious that 

g/3/i j3ane^^. Now expanding the Bose-Einstein function by Taylor series about 
we have the shift of density of the excited bosons 

,\/27rm. 


n 


h 




( 17 ) 


It is seen from equation (17) that the shift of the density vanishes as a —)■ 0 or T —)■ 0, 
which justihes, at least qualitatively, that our derivation is correct. 


D. Casimir amplitude 


For D/\ » 1 the summation over r in equation (6) may be approximated by an 
integration. This approximation leads to 
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where the Casimir amplitude in units of ksT is 
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A 

where Ka denotes the modihed Bessel function. A(T, a) in equation (18) gives a non¬ 
zero magnitude at the critical point and depends only on the boundary condition for a 
specihc dimension. Away from the critical regime it decays exponentially (exponential 
form is implicit in Kd{a)) and thus becomes zero for large D. So, A(T, u) can be 
referred to as the Casimir amplitude in a way analogous to the ideal Bose gas[38]. As 
the IMBG approaches to the point of phase transition i.e. /i —)■ ani^\ a becomes very 



small. In this situation equation (19) reduces to 


= ( 20 ) 

2d--^7r2 / 

Here, we have used the relation Ka{x) ~ (|)" as a; —)■ 0. For d = 3, the Casimir 

amplitude stands as A = —C(3)/(87r) which coincides exactly with the result of the 
ideal Bose gas [38]. 


E. Casimir force 


The Casimir force corresponding to equation (6) is 


f) u„'r (2al) 2 Kd+2(2al) — (2al)2 Kd(2al) 

^ ( 21 ) 

Now for the IBG, a = 0 and then consideration of p ^ 0 in equation (21) gives 
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If we now examine this expression for d = 3, we have 
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(23) 


In order to reexamine the result found in equation (23) and, also to see explicitly 
the behavior of exponential decay we first integrate the third term on the right side 
of equation (6) over r in the three dimensional space. Secondly, we differentiate the 
resulted expression with respect to D. We finally arrive at the following expression for 
the Casimir force 
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where the characteristic decay length 
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In the condensate (T < Tq) fj, = fic = anc- In this case rj ^ oo and equation (24) 
reduces to the form of equation (23), 


ksT C(3) 


Again for T > T^, /i 7 ^ as a result rj is hnite. Then, due to the fact that D jX » 1, 
^ >> 1, so Fc ~ The Casimir force, therefore, shows a power law behavior in 

the condensate, and decays exponentially in the non-condensate (i.e.for T > Tc) with 
a characteristic length rj. It is also noticed from equation ( 21 ) that for T < Tc the 
Casimir force decreases with decreasing value of T and vanishes at T = 0 K. 


F. Critical exponents near the point of Bose condensation 


In order to examine the values of the critical exponent let us start with the charac¬ 
teristic decay length for a IBG for which a = 0. In this case equation (25) yields 

(27) 

The critical exponent is therefore 1/2. This is the same value as shown by others 
[33,39,40]. The critical exponent in the case of d-dimensional space may be obtained 
in the following way. When a —)■ 0, one can easily show for 2 < d < 4 that [1] 

(28) 

where a = /3{/j, — an). A straight forward manipulation leads to (with A'^/ic = a({d) 
and X'^fj, = agd) 
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Now substituting equation (29) into (25) we have the characteristic decay length in the 
d-dimensional conhguration space 

0.5 A 


d pQ d-. 1 


(fl - fic) 


(30) 


So, for the IMBG in d-dimensional space the critical exponent is l/(d-2). For d=3 it 
is just 1 unlike ideal Bose gas for which it is 1/2. The calculated magnitude of the 
exponent agrees well with values evaluated by using a different theoretical method [26]. 
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III. RESULTS AND DISCUSSIONS 


Properties of massive boson gas in the presence of the mean held potential, and con- 
hned between two parallel slabs are presented in this section. This study is performed 
by using the density of states approach with spatial dimension d. Most importantly, 
this approach allows us to study properties of the two phase system, such as the con¬ 
densed and normal phase across the Bose condensation point and, also to examine the 
critical exponents close to the point of phase transition. 

Equation (15) shows the critical temperature, Tc, for the Bose condensation of 
the IMBG. Here Tq is proportional to a factor [ni^^ which is valid for d > 2 
because at d = 2 the Zeta function becomes inhnite and consequently Tc = 0. On 
the other hand, diverges for d = 2. Equation (15), however, demonstrates that the 
critical temperature decreases with increasing d, but for d —)■ cxd the critical temperature 
becomes a constant, Tc = {h/y^TrrmY/ks- 

It is noticed from equation (16) that the fraction of the boson number in the con¬ 
densate varies with the dimensionality in such a way that rio/n increases non-linearly 
with increasing power of d and, as d —)■ cxd all bosons drop to the ground state for 
0 < T < Tc- Equation (16) also shows that at T = Tc no particles are there in 
the ground state and, at T = 0 all bosons fall to it. But in the temperature range 
0 < T < Tc and for hnite d some bosons remain in the excited state. Equation (17) 
shows the shift of density in the lowest order due to the presence of MFP. The sign 
of the shift is negative, that means the density of the excited bosons reduces due to 
the repulsive mean held potential. It is seen from equation (17) that the shift of the 
density vanishes as a —)■ 0, that means if there is no held there is no density shift. It 
is also noticed that as T —)■ 0 shift of density becomes zero. The reason of it is the 
following. At T = 0 K all bosons in the condensed state fall to the ground state as we 
have described before. So no particle remains in the excited state to cause the shift to 
happen. 

Figure 1 shows the Casimir amplitude (in unit of ksT) for jj, = (Eqn. (20)) and 
for jj, ^ He { Eqn. (19)). The former denotes the condensed state and the latter the 
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FIG. 1. d-dependence of the Casimir amplitude, A{d,a)] a is defined in the text. 

normal state. But both properties approach to zero with increasing dimension d. Figure 
illustrates that, in the condensed state (i.e. fx = /ic) the values of A{d, a) is negative and 
large at lower dimension and increases roughly exponentially for the large value of d. 
In the case of non-condensed state the values of the Casimir amplitude are positive and 
large at lower dimension, becomes negative at around d = 3 and, dually approach to 
zero with further increment of d. In the condensed state(/r = /ic ) the Casimir amplitude 
is A(3) = —C(3)/(8 7r), C(3) being the Riemann zeta function. This magnitude is just 
equal to that of the ideal Bose gas [38]. We note here that, besides dimensionality, 
the magnitude of the Casimir amplitude depends on the boundary condition employed. 
For example, A(3) in the case of periodic boundary condition is just eight times larger 
than that of the Neumann or Dirichlet boundary condition [33]. However, equation (20) 
clearly states that A(3) in condensate is not directly affected by the MFP. The reason 
is that, in the condensate [Xc = Cinc which makes z' = 1. Consequently, the Casimir 
amplitude becomes independent of MFP. 

Equation (21) expresses the Casimir force in the case of a IMBG confined between 
two inhnite slabs placed in a d-dimensional conhgurational space. In the limit a —)■ 0 
(i.e. for an ideal Bose gas) this equation reduces to the form of equation (22) which 
hnally takes up the shape of equation (23) for d = 3. This completely agrees with the 
previously established expression for the IBG [14]. This, in turn, justifies that equation 
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(21) correctly represents the Casimir force between two slabs in h-dimensional config¬ 
uration space. The modihed Bessel function Kd in equation (21) decays exponentially, 
consequently the Casimir force decays in the same manner. This decay property can be 
demonstrated more spectacularly by evaluating the Casimir force in a spacial case for 
d = 3 (see equation (24)). Here, equation (24) illustrates that the decay is happening 
with a characteristic decay length r]. We note that the magnitude of rj for a IMBG 
depends on the boundary condition used, for example, = tjd = f]p/‘^ where sub¬ 
script N, D and P denote the Neumann, Dirichlet and periodic boundary conditions, 
respectively [40]. We also note that the decay length is related to the bulk correlation 
length as y/nr] in the case of IMBG and, as —)■ ? 7 p in the case of non-condensed 
phase. This equation also shows that the exponential decay of the Casimir force is also 
related to the separation distance, D, of the two slabs and goes to zero when D is large. 
It is worth mentioning, here, that as /r ^ /ic equation (24) yields the same expression 
for Fc as equation (23). This result also agrees with that of IBG [14]. 

Finally, we turn to the results of critical exponents evaluated near the Bose conden¬ 
sation point. Equation (27) shows that for an ideal Bose gas the critical exponent is 
1/2. On the other hand, equation (30) shows that, in the case of IMBG the exponent 
is l/{d — 2) in the d-dimensional conhguration space which reduces to 1 for d = 3. 
Most importantly, the value of the critical exponent varies as an inverse function of 
dimensionality d. This result agrees with the Endings of others [26,33] although they 
used a difierent theoretical approach. 


IV. CONCLUSION 

The Casimir interaction energy for an IMBG confined between two parallel infinite 
slabs separated by a distance U in a particular direction in a d-dimensional space is 
evaluated, in this report, by using the one particle density of states method. Here the 
Neumann boundary condition is applied, unlike other authors, and the implication of 
other boundary conditions are discussed where appropriate. This study nicely describes 
the amount of density shifts in the lowest order due to the presence of MFP and, to the 
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best of our knowledge this is done for the hrst time. This shift disappears while the 
potential energy goes to zero that means no shift happens in the absensce of a potential. 
The Casimir force obtained by taking derivative of the surface interaction term is found 
to vary as an inverse power law of separation distance D in the condensate and, to decay 
exponentially in the normal state close to the condensation point. The exponent in the 
power law is solely determined by the dimensionality d (~ The Casimir force 

decreases with decreasing value of temperature when T < Tc and vanishes at T = 0. 
On the other hand, in the high temperature limit i.e. for T » Tc the exponential 
factor of the Casimir force approaches to zero. Physically this means that, for T » Tc 
the Bose-Einstein statistics approaches to the Maxwell- Boltzmann one, consequently 
the Casimir interaction term arising due to quantum effect does not appear in the grand 
potential and thus the Casimir force ceases to exist. The critical exponent in the case 
of the IMBG is l/{d — 2) which rednces to 1 in the 3-dimensional conhguration space 
whereas the value of the exponent in the case of IBG is just 1/2. Regarding the critical 
exponent, equation (30) is valid for any value of d except for those of even integers. 
We hnally conclude with the remark that, in case of IMBG, the one particle density 
of states approach is capable to predict the same results as those predicted by the 
rigorous N-body partition fnnction method [33,40] in a rather simpler way. Moreover, 
the present approach correctly describes the amount of density shift which yet to be 
seen from any rigorous approach. 
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